A simple trapped-ion architecture for high-fidelity Toflbli gates 
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We discuss a simple architecture for a quantum TofFoIi gate implemented using three trapped ions. The gate, 
which in principle can be implemented with a single laser-induced operation, is effective under rather general 
conditions and is strikingly robust (within any experimentally realistic range of values) against dephasing, heat- 
ing and random fluctuations of the Hamiltonian parameters. We provide a full characterization of the unitary 
and noise-affected gate using three-qubit quantum process tomography. 

PACS numbers: 03.67 .Lx,32.80.Qk,37.10.Ty 



In the quest for scalability of a quantum computing device, 
the role played by many-qubit gates is quite central. Ade- 
quate sets of gates have been identified that allows for the 
break-down of complex computational networks in simpler 
tasks involving at most two qubits per time 1 1 1 . Unfortu- 
nately, the overhead in terms of the length of corresponding 
quantum circuits (i.e. the number of such elementary opera- 
tions being required) soon overcomes the advantage provided 
by having to manage only two-body interactions. Multi-qubit 
gates exist, able to bypass such a problem by requiring the 
simultaneous conditional evolution of three or more qubits. 
Among these, the Toffoli gate |2| is celebrated for its role in 
phase-estimation and error coiTection protocols |3 1, as well as 
in the quantum factorization algorithm |4|. Remarkably, the 
Toffoli gate has jumped across the field of computing, from 
classical to quantum, playing an important role in schemes 
for reversible classical computation |5|. The recent effort put 
in the task of harnessing a three-qubit Toffoli gate has been 
considerable. On one side, important experimental demon- 
strations on the practical realization of such gate have been 
reported |6, 7|. On the other hand, significant improvement 
in the design of economic ways of implementing an n-qubit 
Toffoli gate have arisen from realizing that less resources are 
needed when using particles that live in higher-dimenional 
Hilbert spaces |8, 9|. Interesting and compact architectures 
using the measurement-based paradigm for information pro- 
cessing have been suggested, in this context fTOl. 

In this paper, we present a protocol for the implementa- 
tion of a three-qubit Toffoli gate in a trapped-ion architec- 
ture. Our scheme exploits an enlarged computational space 
consisting of three-level particles and an ancillary phononic 
mode jn |. While information is encoded only in two elec- 
tronic states of each ion, their third levels are used as conve- 
nient working spaces, similarly to the phononic ancilla. In this 
respect, our algorithm is close in spirit to the work by Ralph et 
al. 1 8 1, although our protocol is different by construction, and 
to further proposals exploiting higher-dimensional particles 
for improved manipulation of the computational states 0. In 
turn, such enlarged computational space allows a consider- 
able reduction in the number of operations used in order to 
implement the three-qubit gate. In fact, our proposal requires 
roughly 44% of the operations needed in the seminal exper- 



imental proof provided by Monz et al. [6]. We thoroughly 
characterize the gate making use of three-qubit quantum pro- 
cess tomography (QPT) |l2l| and reveal its striking robust- 
ness against leakage from the computational space, dephas- 
ing, heating and laser-power fluctuations. We thus provide 
an accurate platform for the experimental implementation of 
such a key gate in the design of scalable computing devices. 

The protocol- Our system consists of A^=3 ions in a lin- 
ear electomagnetic trap. We consider three internal energy 
levels {\lj), \gj), \e j}} in a ladder configuration and the centre- 
of-mass mode a at frequency v. If the ions interact with an 
external classical laser at frequency coi it is possible to couple 
the external and internal degrees of freedom of the ions. In the 
interaction picture with respect to the energy of the internal 
degree of freedom and the vibrational mode, the Hamiltonian 
describing such coupling reads 

i/,(0=(riQ/2)o-L''^^e-'''("''"'''+"*''''''>-'<"''^-'"^''-H/!.c., (1) 

where Q.eR is the Rabi frequency of the transition \a}^\/3} 
(with a,/3-e,g, I), &^^'^^-\a}{f^\, ojap is the corresponding tran- 
sition frequency and rj is the Lamb-Dicke parameter |13|. 
Moreover, we have introduced the phononic annihilation (cre- 
ation) operator a (a^) of the quantized phononic centre-of- 
mass mode. In the following, we will be exploiting the well- 
known flexibility in laser-induced trapped-ion dynamics that 
is achieved by tuning the laser-ion detuning dap-cjJap-^^L fT3l- 
In the Lamb-Dicke regime defined by 77<scl, we can expand 
Eq. ([TJ up to the second order in rj and set dap-v (i.e. we 
tune the laser to the first red sideband of the coupled spin- 
phonon system embodied by a single ion and the centre-of- 
mas mode). This engineers the energy-conserving coupling 
Hl"^\^)-^a&^"^^ -\-h.c., where (-firi{Q./2) and a phononic ex- 
citation is created (destroyed) upon annihilation (creation) of 
a spin quantum. Similarly, the choice 6aj3--v (corresponding 
to the tuning to the first blue sideband) induces the coupling 
H^l^^\i^)-i^a^ &^"^^ +h.c . where spin and phononic excitations 
are simultaneously created or destroyed. We now show how to 
realize a Toffoli gate using the unitary evolution given by the 
Hamiltonian H, , together with properly aiTanged single-qubit 
operations preceding and following the dynamics induced by 
Hr- We will be working in the single-excitation sector of the 
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Hilbert space of the whole ionic string, including the phononic 
mode. Such operations will be required to guarantee that the 
state of the system remains within such subspace. Moreover, 
we want to avoid the presence of correlations between the in- 
ternal degrees of freedom of the string and their vibrational 
one. Therefore, starting from a phononic mode with no exci- 
tations, we have to enforce that when the Toffoli gate is com- 
pleted, the centre-of-mass mode is back to its vacuum state. 

We are now in a position to describe the details of 
our protocol. First, we codify three qubits in the 
internal degrees of freedom of the ions, by using 
the simple encoding scheme (|0i) , |li))=(|^i) , 
(|0,),|l,))=(|g,),|/,» 0=2,3). With this, we con- 
struct an eight-state basis for the three-qubit system as 
S={|000), 1 100), 1010), |1 10), 1001), 1101), -/|01 1), l))i23, 
where we have redefined the last two states so as to include 
an overall phase factor (the choice is made only to simplify 
our calculations). Our protocol begins with the realization 
of the single-qubit operations on the qubit 1 -phononic mode 
system given by (t-7:/20 
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The first transformation in Eq. (j2]) excites the first blue 
sideband for the \g\)^\e\) transition while the remain- 
ing two embody the first red sideband excitations for the 
\ei)*^\h) and |^i)<->|/i) passages. In order to work with the 
same phononic mode, we need laser fields with frequencies 
u'l-LOge=-v,(jUi^-LOie-v,(Ji-ijjgi-v. In the subspace with at 
most a single excitation, it is straightforward to see that the 
composite operation 'R-R^{ji I2(,) R^{k I2() R^{n I2(), oper- 
ated on states having initially no phononic excitations, per- 
forms a logical &^ gate in the space of the phononic mode, 
controlled by the spin state \g\). That is 
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It is important to notice that this operation encodes a logi- 
cal qubit in the single-excitation pair of states {\gi, 1), |ei,0)). 
Therefore, at this stage, one has to consider the attainable 
states of the extended system comprising both the internal and 
external degree of freedom of the ion-string and having a sin- 
gle overall excitation. This is re-interpreted, in our scheme, 
as the computational space of three logical qubits, one of 
which being embodied, in a sort of dual rail encoding, by 
{|OL),|lL))={|^i,l),ki,0)) riTl. We now couple each ion to a 
field of frequency ui-tOeg+v at Rabi frequency Qy. The total 
interaction is described by the Tavis-Cummings model IIT4J 



^TC = Y}-nriajl2ra\ej){gj\ + h.c. 
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We let each of the single-excitation states mentioned above 
evolve under Hie for a time fj- at which we turn Eq. (|4]) off 
and apply the gates forming 'R in reverse order. This sepa- 
rates the state of the centre-of-mass mode from the spin state 



of ion 1. That is, we decode the logical qubit {[Ol),!!/.)) re- 
setting the vibrational mode into |0) and remaining with just 
the ion-string computational states. The vibrational degrees of 
freedom can now be traced out without affecting, in principle, 
the resulting three-qubit gate 

VAf) = ® i23)exp[-(//^)fl-Tcfr](^^ ® 123), (5) 

where ^23 is the identity operators in the tensor-product space 
of qubits 2 and 3. The main idea now is finding an instant 
of time tj such that ^riO is as close as possible to the Tof- 
foH gate f=ii®(i23-|ll)23<ll|)-Ho-;,,i«>|ll)23<ll|. We first 
observe that there are four time-scales associated with the dy- 
namics at hand. Each of them is determined by the inverse of 
the Rabi frequencies of the processes involved in this scheme, 

i.e. 0i23=Mi:,li^^])'^',0i./=Mi^?+i^-)'^' (;=2,3) and 
01=^77^1. The pedices used in such expressions identify the 
qubits that participate to the interaction with the phononic 
mode. It is thus clear that the quest for tj is equivalent to 
the research of a set of suitable single-ion Rabi frequencies 
{Qj) such that our goal is achieved. An extensive numerical 
optimization (performed in the spirit of optimal control the- 
ory) leads to the following relative ratios of coupling strengths 
^1:^2:^^3 = 1: Vl43:16. With this at hand, at the optimal in- 
stant of time given by tT-Tr/rjO.], 'Ut take the following ma- 
trix form [in the ordered computational basis {|/)) .., 8). 
Here ; stands for the decimal-number value of each element 
of the three-ion basis S] 
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Needless to say, other choices can be found for the set of Rabi 
frequencies that achieve a gate close to a Toffoli one. How- 
ever the latter cannot be exactly achieved due to the necessity 
of maximize, simultaneously, eight trigonometric functions of 
incommensurate frequencies. Clearly, the only important pa- 
rameter in our model is the ratio of the Rabi frequencies rather 
than their actual value. The time needed in order to imple- 
ment the whole gate is tc - in I T])[2Yjk=a,b,c^k'^ where 
Q.k is the Rabi frequency of pulse k-a, b, c in Eq. It is 
crucial to stress that in such unitary picture, the implementa- 
tion of Eq. Q for a time tr would be sufficient to implement a 
full Toffoli gate over the logical target qubit (lO^), 1 1/.)), which 
shows the striking economic nature of our proposal. Clearly, 
the use of excited vibrational states would open the protocol 
to the effects of phononic heating and losses. The necessity 
of removing such excitation (so as to make the gate robust) 
motivates the use of the encoding-decoding steps given by K. 

To evaluate the quality of our proposal, the dynamics en- 
compassed by the physical processes described so far should 
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be characterized in a state-independent way. In what follows, 
we take an experiment-inspired approach and use QPT H2J 
as the tool to estimate the performance of the gate. Any com- 
pletely positive A^-qubit map <i>{t) is specified by a set of 4^^ 
orthogonal operators {iC,,,} such that 



(7) 



where p(0) is the initial density matrix of the system and 
we have introduced the A^xA^ time-dependent process ma- 
trix x(t) that incorporates full information on the details of 
the evolution embodied by <!>(?). Pragmatically, this descrip- 
tion is of the utmost usefulness as the map is fully spec- 
ified by considering only a fixed set of operators, whose 
knowledge determines the process matrix. The knowledge 
of the process matrix allows to evaluate the closeness of the 
mechanism under scrutiny to the ideal one corresponding to 
T (with process matrix xt) by means of the gate fidelity 
'^gitG)-Ti[XTX(fc)]- In turn, this is useful to determine the 
average state fidelity associated with gate T Ate)- This pro- 
vides a state-independent indicator of the quality of the gate 
by averaging the fidelity between the ideal and actual out- 
put states over all pure inputs. For a li-dimensional problem 
(^^=2^^) we have [ISJ W s{tcXdTg{tc)+l)l{d+l). InFig.[l]we 
show the representation of the reconstructed process matrix 
in the tensorial operator-basis constructed by considering the 
single-qubit operators {t,&x,-id-y,a-,]. The entries of xitc) 
differ from those of the ideal one by (9(10^^), thus witnessing 
the excellent quality of the gate we have achieved. An addi- 
tional confirmation comes from achieving an average infidelity 
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FIG. 1: (Color online) Reconstructed process matrix corresponding 
to the unitary process 'Uj- The matrix is expressed in the three-qubit 
operator basis formed by the elements {/= 1 , X=(rv, Y= - i&y, Z=&.]. 
We show the moduli of the matrix entries. The differences with re- 
spect to the elements of an ideal gate are (9(10 



i-r.(fG)~io-^ 

Analysis of imperfections.- So far, we have considered only 
unitary evolutions. In order to provide an estimate of the effi- 
ciency of the gate under realistic experimental conditions, we 
should consider the most severe sources of imperfections in 
the ion-trap architecture addressed here [6|. In the following, 
we concentrate on quality-limiting effects of a non-technical 
nature and take into account decoherence of the quantum in- 
formation stored in the phononic mode given by the vibra- 
tional centre-of-mass mode of the ion-string as well as heating 
due to the coupling between the phononic mode and a bath at 
finite temperature. The possibility to tune the Rabi frequency 
and the semiclassical approach used in the ion-pulse dynam- 
ics (implying the use of intense fields in the construction of 'R) 
allow us to neglect the time needed to implement the single- 
qubit operations Eq. Q. We thus consider the master equation 

5,p(f)=-i[^jc,p(f)]-(*f/2)(n-i-l)(a"'"ap(f)+p(0fl'fl-2flp(f)fl') 
- (Kh/2)(aa^ pit) + p(t)aa^ - 2a^pa)—'y[a^a, [a^a,pit)]] 

(8) 

where p(f) is the density matrix of the ionic string including 
the centre of mass mode, k is the heating rate, n is the mean 
number of phononic quanta of the bath at a given tempera- 
ture and y is the dephasing rate. Analogously to the unitary 
case, the dynamical map &h arising from Eq. (|8]l should be 
preceded and followed by the gate. That is, any initial state 
p(0) of the three-ion system evolves until time according to 



p(tc)='R^[&HmpiO)'R^m- 



(9) 



The resulting open-system dynamics implies, in principle, 
leakage from the computational space and thus the spoiling 
of the desired Toffoli gate. In particular, the thermal evolu- 
tion included in our assessment could lead us to abandon the 
subspace where the bosonic mode is in its vacuum state. The 
effective occurrence of such events and their influences over 
the average performance of our operations can be estimated 
by determining again the closeness of the map in Eq. (j9| to 
the ideal transformation T. We have thus used our toolbox for 
QPT to quantify the average gate fidelity for the noise-affected 
evolution and estimate the deviations from ideality. 

We have first considered the effects of heating of the 
phononic mode due to noisy electric potentials on the surface 
of the trap electrodes, resulting in an effective bath at non-zero 
temperature. We have taken (/c"','>'"')=(140, 85)ms and con- 
sidered increasing values of n, taking as reasonable maximum 
number of thermal excitations equal to 5 |6|. In such a worst 
case scenario, we have reconstructed the process matrix x(fc) 
and checked its resemblance to Xt- We have calculated the 
discrepancy \xitc)-XT\ and taken its maximum value per row 
of this matrix (see the main panel of Fig. |2]). The largest devi- 
ation out of the 64 values gathered in this way is a rather small 
^ 2.5x10"^'. In fact, the evaluation of the average gate fidelity 
leads to ^5=0.994855, which is 99.5% the value achieved for 
«=1. The remarkable insensitivity of the scheme to the ef- 
fects of an increased mean phonon number is therefore fully 
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FIG. 2: (Color online) Main panel.- We take the largest entry per row 
in the discrepancy matrix \x(to)-XT\ for "=5, 7/n=10"''. We have 
highlighted the bars corresponding to some of the operator-basis ele- 
ments. Inset.- Lower horizontal axis and circle-shaped points: Aver- 
age state fidelity for the 'lij gate against the dephasing y/Q.. At 7=0 
it is 'Fj=0.999988, while for the larger dephasing rate that we have 
considered we have T's>0.93. Upper horizontal axis and square- 
shaped points: Average gate fidelity for "Wt- against the variance A 
of the uniform distribution determining the amplitudes of laser fluc- 
tuations in our model. The solid lines are only guides to the eye. 



proven. Our analysis thus leads us to conclude that map Sh 
results in an effective dynamics that is approximated in an ex- 
cellent way by p(tc) ~ Pqifc) ® |0)(0|, where Pq(tc) is the 
density matrix of the three-ion system. 

Our next step is the evaluation of the dephasing effects, 
which has been performed by solving the master equation (|8]l 
for n-1, = 140ms and for increasing values of 7. The inset 
of Fig.|2](circle-shaped points) shows the quasi-independence 
of the effective gate from an increase in y by almost one or- 
der of magnitude with respect to the value estimated in re- 
cent cutting-edge experiments (l-;F.5e[10"^, 0.07] for y/Q. e 
[0,12.5]xlO-3) a. 

As remarked above, a key point in our proposal is the main- 
tenance of precise ratios of the Rabi frequencies of the op- 
erations involved in the construction of 'Z/t-. Lasers fluctua- 
tions may, in principle, jeopardize the stability required in the 
scheme we have proposed. Although, experimentally, such 
contributions account for only a negligible percentage of the 
overall gate fidelity [6 |, the central role played by the Rabi 
frequency ratios in our protocol have prompted us to investi- 
gate them thoroughly. We have thus solved Eq. ([8j again, this 
time treating the Q/s as stochastic variables which randomly 
oscillate around the corresponding ideal values. More in de- 
tails, we have taken D.j-Q.j+6Q.j, with 6Q.j/Q.j a uniformly- 
distributed zero-mean variable with variance Ae[l,5]%. Us- 
ing a sample of 500 randomly drawn values of SQ.j and eval- 
uating the corresponding dynamical evolution, we have cal- 
culated the sample-averaged !F.s. In the worst case scenario 
given by A-5% (which grossly overestimates the current ex- 
perimental capabilities), we have achieved an average fidelity 
of ^71% (see the square-shaped points in the inset of Fig.|2]). 

Conclusions.- We have discussed a scheme for the imple- 
mentation of a three-qubit ToffoU gate that requires only 44% 



of the total number of operations needed by a very recent ex- 
perimental demonstration in a trapped-ion system. Despite 
using working space that lies outside the Hilbert space of 
three qubits, our protocol is remarkably robust and affected by 
only negligible leakage from the computational space. More- 
over, it enjoys the gate-catalysis effects provided by the use 
of higher-dimensional information carriers. We have charac- 
terized the performances of the gate from an experimentally 
oriented viewpoint and making use of powerful tools typical 
of QPT. On the one hand, this has allowed an agile quan- 
tification of the overall gate quality. On the other hand, we 
have rigorously estimated the influences of the most relevant 
causes of gate imperfections in the setup at hand, finding quite 
a striking resilience. Economic schemes such as ours will be 
important in the design of experimental-friendly architectures 
for trapped-ion quantum computing that can find a rather non- 
demanding implementation in state-of-the-art settings |i,6. ,16J . 
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